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ABSTRACT 

In order to trace the origin and evolution of carbon in the Galactic disk we have 
determined carbon abundances in 51 nearby F and G dwarf stars. The sample is 
divided into two kinematically distinct subsamples with 35 and 16 stars that are 
representative of the Galactic thin and thick disks, respectively. The analysis is based 
on spectral synthesis of the forbidden [Ci] line at 872.7 nm using spectra of very high 
resolution (R « 220 000) and high signal-to-noise (S/N > 300) that were obtained 
with the CES spectrograph on the ESO 3.6-m telescope on La Silla in Chile. We 
find that [C/Fe] versus [Fe/H] trends for the thin and thick disks are totally merged 
and flat for sub-solar metallicities. The thin disk that extends to higher metallicities 
than the thick disk, shows a shallow decline in [C/Fe] from [Fe/H] w and up to 
[Fe/H] m +0.4. The [C/O] versus [O/H] trends are well separated between the two 
disks (due to differences in the oxygen abundances) and bear a great resemblance to 
the [Fe/O] versus [O/H] trends. Our interpretation of our abundance trends is that 
the sources that are responsible for the carbon enrichment in the Galactic thin and 
thick disks have operated on a time-scale very similar to those that are responsible 
for the Fe and Y enrichment (i.e., SNIa and AGB stars, respectively). We further 
note that there exist other observational data in the literature that favour massive 
stars as the main sources for carbon. In order to match our carbon trends, we believe 
that the carbon yields from massive stars then must be very dependent on metallicity 
for the C, Fe, and Y trends to be so finely tuned in the two disk populations. Such 
metallicity dependent yields are no longer supported by the new stellar models in 
the recent literature. For the Galaxy we hence conclude that the carbon enrichment 
at metallicities typical of the disk is mainly due to low and intermediate mass stars, 
while massive stars are still the main carbon contributor at low metallicities (halo and 
metal-poor thick disk). 

Key words: stars: abundances - stars: kinematics - Galaxy: abundances - Galaxy: 
disc - Galaxy: evolution 



1 INTRODUCTION 

Next to hydrogen, helium and oxygen, carbon is the most 
abundant element in the Universe and plays an important 
role in the chemical evolution of galaxies. Although the nu- 
clear process (helium burning) that is responsible for the 
for mation of carbon atoms is well-know n (see, e.g., reviews 
bv lWallerstein et alJll997l:lArnettll2004 . it is unclear which 
objects that contribute to the carbon enrichment of the in- 
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terstellar medium. This has been m uch debated during re- 
cent years; iGustafsson et ail il999T) conclude that carbon is 
mainly contributed f rom massive, metal-rich stars, and not 
from low- mass stars; IChiappini et al ] J2003allbll find strong 
indications that carbon is p roduced in low and intermediate 
mass stars; IShi et all i2002tl find that carbon is contributed 
by super-winds from massive stars in the early stages of 
disk formation in the Galaxy, while at later stages signifi- 
cant amounts of carbon are contributed by low-mass stars. 
Other studies that favour a mixture of low mass, interme - 
diate mass, and high mass stars include l Liang e t al.| l|200l[); 
iGavilan et all (120051) : ICarigi et al J (120051) while iHenrv et all 
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(2000) and lAkerman et alJ i2004T) suggest massive stars to 
be the main source. 

Only few previous stud ies lAjioLers^pji_fcEc lvardssonl 
119941: IGustafsson et alJll999t iTakeda fc Hondal 120051) have 
based their carbon abundances on the forbidden [C i] line 
at 872.7 nm. This line is highly insen sitive to errors in the 
stellar model atmosphere (see, e.g., lAsplund et alJ 120051) 
and gives carbon abundances of high accuracy in solar- 
type stars, even when using one-dimensional (1-D) models 
under the assumption of local thermodynamic equilibrium 
(LTE). Unfortunately, the line is very weak and is blended 
with a weak feature tha t is likely to be an Fel line (e.g. 
lLambert fc Swings! Il967l) . So, in order to achieve accurate 
abundances from this line it is desirable to have spectra of 
high- resolution and high signal-to-n o ise (S /N) and it is es- 
sentially only the IGustafsson et all (I1999T) study that has 
reliable carbon trends based on the forbidden line (see dis- 
cussion in Sect. 14.51 . They did, however, not include the 
blending Fe I line in their analysis. 

As a step to achieve a deeper understanding of the sites 
of carbon formation and the chemical histories of the Galac- 
tic disks we have obtained high-resolution, high S/N spectra 
of the forbidden [C i] line at 872.7 nm for 51 dwarf stars 
in the solar neighbourhood. These stars are already well 
studied, have well determin ed oxygen and iron abundances , 
and have known kinematics feensbv et alj2003ll20oll2005ll . 
That the kinematics are known is important as this makes 
it possible to properly investigate the carbon trends in the 
two disks, independent of each other, and to perform a dif- 
ferential analysis of the carbon trends in the thin and thick 
disks in order to investigate the differences and similarities in 
the disks' chemical histories. The latter investigation could, 
combined with other evidence, be important for understand- 
ing the time-scales involved in the chemical enrichment of 
the gas that today make up what we observe as the stars in 
the Galactic thick disk. 

In Sect. |21 we shortly describe the stellar sample and 
the observations. Sect. |3] discusses the data reduction and 
associated problems and Sect. [3] discusses the abundance 
analysis and estimate random and systematic errors in the 
derived carbon abundances. In Sect. [K] we describe our re- 
sulting carbon trends and we use our results to trace the 
possible sources of carbon as well as discuss the chemical 
evolution of the Galactic disks. Finally, Sect. |S| concludes 
with a summary of our findings. 



2 STELLAR SAMPLE AND OBSERVATIONS 

The stellar sample (listed in Table contains 51 F and 
G dwarf stars (35 thin disk and 16 thick disk stars) which 
form a subset of the samples in our previous s tudies (in total 
102 stars, see lBensbv et all200ll2004l2005h . The selection 
criteria used to assign a star to either the thin or the thick 
disk as well as further details and discussions about the kine- 
matic properties of the st ars in the sample can be found in 
iBensbv et all fcOOl l2005l) and in Sect.ETTl 

Observations were carried out by TB during six nights 
in September 2004 with the Coude Echelle Spectrograph 
(CES) on the ESO 3.6-m telescope on La Silla in Chile. 
All spectra were obtained with the highest resolution R ~ 
220 000. Although our exposure times were calculated to 




Figure 1. Top: A normalised flatficld frame. Middle: A "raw" 
science frame (in this case a solar spectrum) Lower: The flat- 
fielded science frame. The strong Si I line at 872.8 nm, in whose 
left wing the [Cl] line is located, has been marked. The full 
CCD frames are 4096x1024 pixels. Here we show the central 
parts (approximately 2000x500 pixels in size) from ~ 871.6 nm 
to ~ 874.0 nm. 
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Figure 2. The average profile (based on 4096 rows) perpendicular 
to the dispersion direction. Note that the background level is 
higher on the right hand side (pixel value ~ 240) than on the left 
hand side (pixel value ~ 220). 



give spectra with signal-to-noise ratios that exceed at least 
350 this was unfortunately rarely the case in the reduced 
spectra. The reason for this degradation in quality of the 
reduced spectra is due to unexpectedly strong fringing pat- 
terns from/in the CCD detector in this wavelength region. 
The treatment of this feature is discussed in more detail in 
Sect. El 

Further, long exposures were split into shorter expo- 
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Figure 3. This comparison shows the efficiency of the Hatfield 
in removing the fringing pattern from a "raw" image. The 1- 
D spectra shown here are the same as those in the two bottom 
panels in Fig.fflwith the exception that the reduced and flatfielded 
spectrum has been filtered for cosmic ray hits (which was not the 
case in Fig.0. The "raw" spectrum has also been normalised and 
vertically shifted for this plot. 



sures of maximum 30 minutes. Solar spectra were obtained 
each day by observing the day-sky approximately one hour 
before sunset. Flatfield and bias exposures were taken at the 
beginning and at the end of the nights and additional flat- 
fields at least once more during the nights. Wavelength cali- 
bration spectra were obtained from a thorium-argon hollow 
cathode lamp. The last four nights we also included expo- 
sures of rapidly rotating B stars that were used in the reduc- 
tion procedure to minimise the residuals from the fringing 
pattern. 



3 DATA REDUCTION 

The top panel in Fig. shows a part of the CCD frame 
for a (normalised) flatfield exposure and the middle panel a 
solar spectrum. As can be seen they are severely affected by 
fringing. However, this fringing pattern was stable and did 
not change over the six nights the observations lasted. As 
the removal of this pattern is extremely important for the 
quality of our results we give a more detailed description of 
the reduction process. 

The spectra were reduced using standard routines in 
the MIDAS 1 software package and consisted of the follow- 
ing steps: 

1) One master bias frame and several (typically 4-5) master 
flatfield frames were constructed for each of the six nights 



ESO-MIDAS is the acronym for the European Southern Obser- 
vatory Munich Image Data Analysis System which is developed 
and maintained by the European Southern Observatory. 



Figure 4. The top spectrum shows the B star template. It has 
been shifted vertically for this plot. Also shown is a reduced solar 
spectrum (dotted line, same spectrum as in Fig.[3J and the same 
solar spectrum but divided by the normalised B star template 
(solid line). 



by averaging multiple exposures. Bias subtraction for the 
flatfields were done after the bias value had been corrected 
for the difference between its overscan value and that of the 
flatfield frame. Since the main method for removing fringing 
patterns is through the flatfielding process care was taken 
when normalising the flatfield frames. As proposed to be the 
best method in the CES manual 2 we normalised the flat- 
field row by row, i.e., a low-order polynomial was fitted to 
each row (i.e., along the dispersion direction) in the flatfield 
frame. This process is necessary since the flatfields are ob- 
tained through the fibre and the image slicer which results 
in a flatfield profile that is composed of the same 12 slices as 
the science exposures (see Fig. A plain division by this 
"non-flat" flatfield would result in too high weights given to 
the low-flux minima regions between the slices (i.e. dividing 
low numbers by low numbers). A normalised flatfield frame 
can be seen in the top panel in Fig. Q 

2) The science frames were first corrected for the bias level. 
Before dividing by a flatfield frame we also corrected the 
frames for the fact that the background level is higher on 
the right hand side than on the left hand side of the profile 
perpendicular to the dispersion direction (see Fig.^l. This is 
due to scattered light originating within the CES and do not 
contribute to the fringing pattern. If it is not corrected for 
before the flatfielding the division by the flatfield will over- 
compensate the pattern. We therefore fitted linear functions 
between the left hand side (approximate pixel number 400) 
and the right hand side (approximate pixel number 800) and 
made a subtraction so that the two sides were level. This 
correction had a very positive effect on the quality of the re- 



Available at http : //www. Is . eso . org/lasilla/sciops/3p6/ces/ 
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Table 1. Stellar names, chemical abundances, and stellar parameters. Abundances are given relative to our solar abundances. [C/H] 
abundances derived without the blending Fel line are given in italics in parenthesis. Metallicities ([Fe/H]), effective temperatures (T c g), 
surface gravities (logg), microturbulences (£t), and space velocities (J/lsr? ^LSRi Wlsr) are taken from Bensby et al. (2003, 2005). 
Oxygen abundances ([O/H]) are taken from Bensby et al. (2004b) and are generally based on the forbidden [Ol] line at 630. Onm unless 
otherwise indicated ("*" means that the permitted Ol triplet lines at A sa 777 nm were used and "**" means that the forbidden [Ol] 
line at A = 636.3 nm was used). fR_x is the macroturbulent broadening (RAD-TAN) that we used (these values have no "real physical" 
meaning and are only listed for completeness). 



HIP 


HD 




\f /HI 




[O/H] 


[C/O] 


Teff 


l°g9 


St 


^R— T 


t-'LSR 


Vlsr 


Wlsr. 
















[K] 


[cgs] 


[km s" 1 ] 


[km s^ 1 ] 


[kms -1 ] 


[kins -1 ] 


[kms" 1 ] 












THIN DISK STARS 














910 


693 


-0.36 


-0.33 


(-0.34) 






6220 


4.07 


1.43 


5.09 


29.0 


-7.8 


-11.1 


2235 


2454 


-0.28 


-0.11 


(-0.13) 






6645 


4.17 


1.75 


8.33 


23.2 


-26.2 


-6.7 


2787 


3229 


-0.11 


-0.05 


(-0.08) 






6620 


3.86 


1.70 


5.44 


-16.5 


-21.7 


-12.9 


3142 


3735 


-0.45 


-0.37 


(-0.39) 


-0.25 


-0.12 


6100 


4.07 


1.50 


4.13 


-35.5 


-45.8 


19.9 


3909 


4813 


-0.06 


-0.02 


(-0.03) 


-0.11 


0.09 


6270 


4.41 


1.12 


3.85 


31.4 


2.5 


-4.8 


5862 


7570 


0.17 


0.21 


(0.20) 


0.06 


0.15 


6100 


4.26 


1.10 


4.96 


-34.1 


-17.1 


-6.0 


7276 


9562 


0.20 


0.19 


(0.18) 


0.14 


0.05 


5930 


3.99 


1.35 


4.38 


1.1 


-21.1 


19.5 


9085 


12042 


-0.31 


-0.20 


(-0.22) 


*-0.12 


-0.08 


6200 


4.25 


1.30 


5.43 


-42.5 


-5.1 


1.0 


10306 


13555 


-0.17 


-0.03 


(-0.05) 






6560 


4.04 


1.75 


8.20 


-10.3 


-6.8 


11.2 


12611 


17006 


0.26 


0.18 


(0.22) 


**0.14 


0.04 


5250 


3.66 


1.35 


3.47 


22.0 


-13.5 


0.1 


12653 


17051 


0.14 


0.15 


(0.14) 


0.04 


0.11 


6150 


4.37 


1.25 


6.67 


-21.3 


-12.8 


-3.1 


14954 


19994 


0.19 


0.21 


(0.20) 


0.11 


0.10 


6240 


4.10 


1.60 


9.24 


-9.7 


-14.2 


1.2 


15131 


20407 


-0.52 


-0.39 


(-0.40) 


-0.26 


-0.13 


5834 


4.35 


1.00 


2.97 


3.6 


20.2 


-9.6 


17378 


23249 


0.24 


0.10 


(0.18) 


0.13 


-0.03 


5020 


3.73 


0.80 


3.15 


-4.8 


31.6 


18.5 


23941 


33256 


-0.30 


-0.20 


(-0.22) 


*-0.23 


0.03 


6427 


4.04 


1.90 


10.0 


-2.0 


-1.9 


8.1 


24829 


35072 


0.06 


0.09 


(0.07) 


0.03 


0.06 


6360 


3.93 


1.70 


4.45 


-36.8 


-26.3 


-15.5 


29271 


43834 


0.10 


0.08 


(0.11) 


0.02 


0.06 


5550 


4.38 


0.80 


2.63 


29.5 


-27.2 


-5.8 


78955 


144585 


0.33 


0.21 


(0.22) 


0.17 


0.04 


5880 


4.22 


1.12 


3.16 


-16.6 


-16.2 


25.4 


80337 


147513 


0.03 


-0.09 


(-0.07) 


0.03 


-0.12 


5880 


4.49 


1.10 


3.08 


24.0 


4.0 


5.8 


80686 


147584 


-0.06 


-0.13 


(-0.13) 


-0.05 


-0.08 


6090 


4.45 


1.01 


3.99 


18.5 


14.3 


2.5 


84551 


156098 


0.12 


0.07 


(0.05) 


0.12 


-0.05 


6475 


3.79 


2.00 


7.61 


-31.9 


-14.7 


16.8 


84636 


156365 


0.23 


0.21 


(0.21) 


0.15 


0.06 


5820 


3.91 


1.30 


4.18 


1.7 


3.8 


-21.4 


86796 


160691 


0.32 


0.25 


(0.26) 


0.14 


0.11 


5800 


4.30 


1.05 


3.44 


-5.0 


-2.7 


3.5 


87523 


162396 


-0.40 


-0.33 


(-0.35) 


-0.31 


-0.02 


6070 


4.07 


1.36 


3.72 


-14.9 


-5.2 


-23.7 


90485 


169830 


0.12 


0.18 


(0.17) 


**0.06 


0.12 


6339 


4.05 


1.55 


4.42 


-6.1 


6.5 


10.9 


91438 


172051 


-0.24 


-0.30 


(-0.28) 


-0.18 


-0.12 


5580 


4.42 


0.55 


2.27 


47.1 


3.0 


2.9 


99240 


190248 


0.37 


0.25 


(0.29) 


0.19 


0.24 


5585 


4.26 


0.98 


3.08 


-38.2 


-8.5 


-8.2 


103682 


199960 


0.27 


0.25 


(0.25) 


0.14 


0.11 


5940 


4.26 


1.25 


3.70 


3.8 


-18.1 


3.6 


105858 


203608 


-0.73 


-0.60 


(-0.62) 


-0.46 


-0.14 


6067 


4.27 


1.17 


3.50 


-3.8 


49.5 


13.7 


109378 


210277 


0.22 


0.10 


(0.15) 


0.16 


-0.06 


5500 


4.30 


0.78 


2.60 


13.6 


-45.5 


1.8 


110341 


211976 


-0.17 


-0.14 


(-0.15) 


-0.11 


-0.03 


6500 


4.29 


1.70 


6.64 


4.7 


11.5 


0.7 


113137 


216437 


0.22 


0.22 


(0.22) 


0.13 


-0.09 


5800 


4.10 


1.16 


3.99 


14.1 


15.0 


5.2 


113357 


217014 


0.20 


0.13 


(0.15) 


0.10 


0.03 


5789 


4.34 


1.00 


3.23 


-5.2 


-23.2 


22.0 


113421 


217107 


0.35 


0.27 


(0.31) 


0.16 


0.11 


5620 


4.29 


0.97 


2.75 


8.3 


-3.8 


18.3 


117880 


224022 


0.12 


0.19 


(0.17) 


0.05 


0.14 


6100 


4.21 


1.30 


5.36 


-36.2 


-11.4 


2.8 












THICK DISK STARS 












3086 


3628 


-0.11 


0.04 


(0.03) 


0.09 


-0.05 


5840 


4.15 


1.15 


2.99 


-159.4 


-48.4 


53.6 


3497 


4308 


-0.33 


-0.25 


(-0.25) 


-0.05 


-0.20 


5636 


4.30 


0.80 


2.97 


60.0 


-103.7 


-19.0 


5315 


6734 


-0.42 


-0.42 


(-0.38) 


0.00 


-0.42 


5030 


3.46 


0.86 


2.41 


60.3 


-118.0 


46.3 


14086 


18907 


-0.59 


-0.53 


(-0.52) 


-0.16 


-0.37 


5110 


3.51 


0.87 


2.54 


18.6 


-78.1 


-13.1 


15510 


20794 


-0.41 


-0.38 


(-0.35) 


-0.01 


-0.37 


5480 


4.43 


0.75 


2.04 


-69.4 


-89.3 


-24.0 


17147 


22879 


-0.84 


-0.76 


(-0.77) 


-0.32 


-0.44 


5920 


4.33 


1.20 


2.29 


-99.0 


-80.4 


-37.0 


18235 


24616 


-0.71 


-0.72 


(-0.71) 


-0.24 


-0.48 


5000 


3.13 


0.95 


2.97 


-17.0 


-158.1 


-20.1 


75181 


136352 


-0.34 


-0.20 


(-0.21) 


0.00 


-0.20 


5650 


4.30 


0.78 


2.58 


-109.5 


-41.4 


43.3 


79137 


145148 


0.30 


0.12 


(0.24) 


0.20 


-0.08 


4900 


3.62 


0.60 


3.03 


81.5 


-49.8 


-62.6 


82588 


152391 


-0.02 


-0.11 


(-0.02) 


0.01 


-0.12 


5470 


4.55 


0.90 


3.32 


94.1 


-105.5 


16.6 


88622 


165401 


-0.46 


-0.44 


(-0.44) 


-0.05 


-0.39 


5720 


4.35 


0.80 


3.17 


-68.6 


-84.7 


-31.7 


96124 


183877 


-0.20 


-0.07 


(-0.06) 


0.08 


-0.15 


5590 


4.37 


0.78 


2.27 


-28.0 


-87.4 


-15.4 


103458 


199288 


-0.65 


-0.56 


(-0.57) 


-0.29 


-0.27 


5780 


4.30 


0.90 


2.60 


31.3 


-96.3 


52.9 


108736 


208998 


-0.38 


-0.30 


(-0.31) 


-0.01 


-0.29 


5890 


4.24 


1.05 


3.00 


-15.6 


-72.9 


44.6 


109821 


210918 


-0.08 


-0.08 


(-0.08) 


0.10 


-0.18 


5800 


4.29 


1.05 


2.72 


-36.9 


-86.9 


-2.5 


118115 


224383 


-0.01 


0.07 


(0.07) 


0.04 


0.03 


5800 


4.30 


1.00 


3.05 


-65.0 


-77.9 


3.3 



© 0000 RAS, MNRAS 00Q.HHl31 



The origin and chemical evolution of carbon in the Galactic thin and thick disks 5 



duced spectra for stars with long exposure times as the dif- 
ference between the two sides relative to the total flux of the 
object seem to increase with exposure time. For the flatfield 
exposures the difference between left and right hand sides 
is practically zero. After this correction we divided the sci- 
ence frames with the normalised flatfield frames. The middle 
panel in Fig.^shows a "raw" science frame and the bottom 
panel in Fig. an example of a frame that has been cor- 
rected for the straylight pedestal and then divided by the 
normalised flatfield frame shown in the top panel in Fig. 
In this way we were able to quite efficiently remove much 
of the fringing pattern which is further illustrated in Fig. 
The final spectra were also filtered for cosmic ray hits (using 
the MIDAS task filter/cosmic) before dividing with the 
flatfield. 

3) One-dimensional spectra were obtained by summing the 
fluxes of all 12 slices (including the regions in between the 
slices, see Fig. The spectra were then wavelength cali- 
brated and re-binned to constant steps in wavelength. Fig- 
ure [3 shows the reduced solar spectrum. 

4) We also observed several rapidly rotating B stars. Since 
their spectral lines are all smeared out their spectra should 
more or less form a featureless continuum. However, the re- 
duced spectra for our 19 B stars have a structure that show 
very little variance when compared to each other (apart from 
variations in the signal-to-noise ratio). This makes it most 
likely that this structure is the remains of an "imperfect" re- 
moval of the fringing and what we see in the B star spectra 
should without any doubt also be present in all our other 
reduced spectra. In order to further reduce the residuals 
from the fringing we therefore divided all spectra with a nor- 
malised B star template. This B star template was derived 
as a weighted average of all 19 B star spectra. Figure^Jshows 
the normalised B star template and the solar spectrum be- 
fore and after division with it. 

5) For stars that have multiple exposures we then co-added 
the individual spectra to one spectrum. Finally, the spectra 
were normalised by fitting a low-order Legendre polynomial 
to the continuum using the IRAF 3 task CONTINUUM. 



4 ABUNDANCE DETERMINATION 
4.1 Stellar atmospheres 

We use one-dimensional, plane-parallel LTE stellar model 
atmospheres that were calculated with the Uppsala 
marcs code llGustafsson et al.ll975l : lEdvardsson et alll99^: 
lAsplund et allll997l) . Surface gravities (logg) were deter- 
mined from Hipparcos parallaxes, metallicities ([Fe/H]) were 
determined using Fel lines, effective temperatures (T e g) 
were determined by requiring that the abundances derived 
from Fel lines with different excitation energies should all 
yield the same [Fe/H], and the microturbulence parameter 
(£t) by requiring all Fel lines should yield the same abun- 
dances independent of line strength (logW\/ X). All these 



3 IRAF is distributed by National Optical Astronomy Observa- 
tories, operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foun- 
dation, USA. 



Table 2. Atomic data. The columns give: element and ioni- 
sation stage, wavelength, lower excitation energy, log gf- value, 
method for collisional broadening (see Sect. |4~2"1 . correction fac- 
tor to the classical Unsold broadening, damping constants, and 
references for the log gf- values: H93 = Hibbert et al. (2003); G97 
= Galavis et al. (1997); K93 = Kurucz (1993); BFL03 = Bensby 
et al. (2003). (Note: see discussion in Sect. 4.4 and in Allende Pri- 
eto et al. (2002) regarding the Fe I line and why log gf = —3.93 is 
likely an overestimate). The full line list that was used in the cal- 
culations of the synthetic spectra can be reproduced by querying 
the VALD database Ihttp: //www. as tro .uu. se/~vald) for lines 
in the region 872.0 -873.5 nm. 
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parameters (listed in Table [D were de termined in our previ- 
ous studies (Bcnsbv ct al. 2003, 2005) wherein the iterative 
process to tune the stellar parameters also is fully described. 

The [Fe/H] values listed in Tableware given relative to 
the solar Fe abundance that we derived using solar spec- 
tra obtaine d with the same spectro graphs as the stellar 
spectra (see iBensbv et al J 120051 120051) . In these studies we 
also compare our derived [Fe/H] values as derived from Fel 
lines to the [Fe/H] values derived from Fell lines (these we 
did not use in the tuning of the stellar paramet ers). Feu 
lines are more robust indicators than Fel lines fsee lAsplundl 
120051) since they are normally not affected by departures 
from LTE. As we generally find good agreement between the 
[Fe/H] value s derived from Fel and Fell lines ; respectively, 
(see F ig. 7 in lBensbv et aljl2003l and Fig. 2 in lBensbv et alJ 
2005) this indicates that for the range of parameters that 
the stars in our investigations span the effects of departures 
from LTE are small, if at all measurable. 



4.2 Synthetic spectra, atomic data, and line 
broadening 

The synthetic spectra were calculated with the Uppsala 
spectrum software. As input it needs a stellar atmosphere 
model and a list of spectral lines with atomic data. To 
properly reproduce t he region around the [C i] line we used 
the VALD database Ku pka et alJll99*S IRvabchikova et al] 
119991 iPiskunov et alJll99Sft to extract all lines in the re- 
gion 872.0 nm to 873.5 nm. Table [5] lists the atomic data 
for the [Ci] line at 872.7126 nm (log gf = -8.136, theo- 
retical from iHibbert et"aT1 Il993l : iGalavis et aD Il997l) . the 
blending Fel line at 872 .7132 nm (log gf = -3.93, theoret- 
ical from lKuruc"3ll993l) . and the Si I line at 872.8010 nm 
(log gf = —0.42, astrophysical from IfBensbv et alJ 120031) 
whose left wing form the continuum level for the [C 1] line. 

The broadening of atomic lines by radiation damp- 
ing was considered in the determination of the abundances 
and the damping constants (7 ra d) for the different lines 
were taken from the VALD database. Collisional broaden- 
ing, or Van der Waals broadening, by hydrogen atoms was 
also considered. The widt h cross-section for the Si I line at 
872.8010 nm is taken from lBarklem et alJ JT998) (indicated 
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872.74 872.76 872.78 872.8 872.82 872.84 
Wavelength [nm] 

Figure 5. Determination of the broadening due to macrotur- 
bulence. (see Sect. 4.2). This plot shows the solar spectrum for 
which we found a best matching RAD-TAN profile with a width 
of 3.11 km s -1 . The synthetic spectra shown in the two bottom 
panels have RAD-TAN widths from 1.5 to 5.0kms _1 . 




872.66 872.68 872.7 872.72 872.74 872.76 
Wavelength [nm] 



Figure 6. The [C i] line in the solar spectrum. The dots are 
the observed spectrum and the thick solid line the best fit repre- 
senting a solar carbon abundance of logeQ(C) = 8.41. The thin 
solid lines represents 5 different carbon abundances from 8.34 to 
8.46 dex in steps of 0.03 dex. The dashed line shows the contribu- 
tion from the blending Fe I line. The wavelength positions of the 
[C i] and the Fe I lines are also indicated. 



by a "C" in Table |5J and for the other two lines we apply 
the correction term (#76 ) to the classical Unsold approxima- 
tion for the Van der Waals damping (indicated by a "U" in 
Table HJ. 

The synthetic spectrum is then convolved with a line 
profile to reproduce the instrumental broadening and then 
with a another line profile to reproduce the broadening due 
to large-scale motions in the stellar atmosphere. The width 
of the instrument broadening profile is set by the resolution 
of the spectra and was in our case 0.0042 nm. This value was 
determined by measuring the FWHMs in the ThAr spec- 
tra and corresponds to a spectral resolving power of about 
R ~ 208 000 at 873.0 nm (R = A/ A A = 873.0/0.0042). 
The large-scale motions include macroturbulence (that has 
a radial-tangential, "RAD-TAN" , profile) as well as the line- 
of-sight component (v-sini) of the stellar rotational velocity. 
The contribution from v ■ sin i is usually small in our types of 
stars (F and G d warf stars) in comparison to the macrotur- 
bulence (see, e.g.. lGravlll992h . To determine the widths of 
the RAD-TAN profiles we used the quite strong Si I line at 
872.8 nm. First we usually had to change the log gf- value of 
the Si I line in order to get the line strength correct. We then 
made 20 different synthetic spectra with different RAD-TAN 
broadenings from to 10 kms -1 . The best fit was found by 
minimising an un-normalised x 2 -function. As we here only 
want to find the the best solution with a given set of free 
parameters we need not normalise our x 2 s- An example of 
this is shown in Fig. 



4.3 Carbon abundances 

The forbidden carbon line is located in the left wing of 
a Si I line a nd is blended by a line that most likely is 
an Fe I line jLambert fe Swings! Il967l) . The contribution 
from this blending line to the joint [C i]-Fe I line profile 
is generally negligible at sub-solar metallicities while it in 
the Sun is estimated to b e between 0.1pm and 0.5 pm 
jAllende Prieto et al.lf2002l) . While most previous studies 
have chosen to neglect this line in their abundance analyses 
due to its assumed weakness and uncertain atomic data we 
have included it as our sample contains many stars at super- 
solar [Fe/H] where such a blend will grow in size. Since thi s 
Fel line only has a calculated log gf- value jKurucjll99l) . 
and could therefore be associated with large errors, we must 
be careful and make extensive checks. For comparison pur- 
poses we will therefore also determine carbon abundances 
neglecting the Fel blend. 

By producing a set of synthetic spectra with carbon 
abundances varying in steps of 0.03 dex we then minimised 
an un-normalised x 2 -function to find the best fitting syn- 
thetic spectrum and hence the carbon abundance (in the 
same way as we determined the macroturbulence broaden- 
ing, see Fig. 0. For seven stars (HIP 10798, HIP 72673, 
HIP 81520, HIP 102264, HIP 107975, HIP 3185, and 
HIP 109450) we did not determine carbon abundances due 
to that their spectra showed defects probably as a result of 
an insufficient removal of the fringing pattern. 

Figure |S| shows an example of the observed and the 
synthetic spectra for the Sun. From our analysis of the 
solar spectrum (using a MARCS model atmosphere with 
T cff = 5777 K, logg = 4.44, & = 0.85 kms" 1 , and RAD- 
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Figure 7. The effect on [C/H] when not including the Fel blend when synthesising the [C i] line. The difference between the two 
derivations of carbon are shown in a) versus effective temperature; b) versus surface gravity; and c) versus metallicity. d) then shows 
the [C/Fe] versus [Fe/H] trend for all stars; e) [C/Fe] versus [Fe/H] for all stars when neglecting the Fel blend; and f) [C/Fe] versus 
[Fe/H] (when including the Fel blend) for stars that have T e ff > 5700 K. 



TAN profile of 3.11 km s -1 , see Fig |SJ we get a solar car- 
bon abundance of log£0(C)8.41 when the blending Fel line 
is taken into account and log e©(C) = 8.44 when it is ne- 
glecte d. Both are in good agreem ent with the recent analy- 
sis bv lAllende Prieto et all (120021) and lAsolund et al.l feOOBh 
who found a best value of log 60(C) = 8.39 using 3-D stel- 
lar atmosphere models and loge0(C) = 8.41 when using 
a MARCS model atmosphere. These studies did, however, 
not include the blending Fel line in their calculations. The 
reason for that we see a difference of 0.03 d ex (8.44 -8.41) 
between our and lAllende Prieto et all i2002T) analysis, both 
using MARCS atmosphere models, is probably due to un- 
certainties in the placement of the continuum and/or dif- 
ferences between the spectra (that could originate from the 
data reduction process, see Sect. El . 

Our final [C/H] abundances, both with and without 
the Fel blend taken into account, are listed in Table 
(given relative to our solar abundance log e©(C) = 8.41 and 
loge©(C) = 8.44, respectively). 



4.4 Fel blend or no Fel blend? 

How does the inclusion of the Fel blend influence our fi- 
nal carbon abundances and carbon trends? In Figs. EK-c we 
show the difference in [C/H] versus Teff, logg, and [Fe/H], 
respectively. There is a clear trend with T e fs; the lower the 
temperature the larger is the effect on [C/H] when neglecting 
the blend. That the difference in [C/H] generally is negative 
below 5700 K and positive at higher temperatures is due to 
that our abundance analysis is differential to the Sun. At 
effective temperatures higher than ~ 5700 K the difference 
in [C/H] is always smaller than ~ 0.03 dex no matter if the 



stars have high or low metallicities. No absolute trend with 
[Fc/H] can be seen (there is still a number of stars at high 
metallicities where the difference is negligible). From this we 
conclude that it is for a combination of low T e g (and often 
low logg), and high [Fe/H] that we need to be extra careful 
and correctly include the Fel blend. 

In Figs. [TJi-e we show the effect of the Fe I blend on the 
final [C/Fe] versus [Fe/H] trend. Apart from a somewhat 
lower spread in [C/Fe] when neglecting the blend the trends 
are essentially identical: both have a flat [C/Fe] trend at 
[Fe/H] < and at [Fe/H] > the [C/Fe] ratio declines when 
going to higher metallicities. Figure |7f shows the trend (de- 
rived with the blend) but for stars that have T e s > 5700 K 
and for which the effect on [C/Fe] due to an erroneous treat- 
ment of the blend should be minimal (see 

Fig-El)- The trend 

seen in Figs. |7Jl-e still persist. 

Summarising, we conclude that it is important to in- 
clude the Fe I blend for stars that have low effective temper- 
atures (below 5700 K) and then especially at high [Fe/H]. 
It is, however, important to remember that the log q /-value 
that is available for the Fe I line is a theoretical one <lKurucz| 
1993) and it is possible that the calculated line strength 
for the blend is either over- or underestimated. It would 
therefore be very valuable to have an accurate experimental 
log g /-value for this Fel line so that it is possible to prop- 
erly calculate its c ontribution to the joint [C l]-Fe I line pro- 
file. We note that lLambert fe RieJ <ll977f) tried to estimate 
the log gf- value for this line by ensuring that the same C 
abundance is derived from the [C 1] line as from C2 and CH 
lines. They arrived at a value of log gf w —3.6 which would 
lead to a higher Fe I contribution to the joint [C 1] - Fe I line 
profile at 872.7 nm than indicated by log gf = —3.93 from 
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Figure 8. [C/H] versus [Fe/H] for our and the stellar samples 
from Gustafsson et al. (1999), Andersson & Edvardsson (1994), 
and Takeda & Honda (2005). Linear regression lines are also plot- 
ted, see Eqs. For our stellar sample we also plot the abun- 
dances that were derived when neglecting the Fel blend (open 
circles and dotted line). 

iKurucd 1^)93). By looking at other Fel l ines b elonging to 
the same multiplet lAllende Prieto et alJ (J20O2) concluded 
in their analysis of the solar spectrum that the contribution 
should be significantly lower and hence that log gf — —3.93 
would overestimate the line strength of the blending Fe I line. 
Allende Prieto et al.'s (2002) conclusion means that our car- 
bon abundances might be somewhat undere stimated at hig h 
[Fe/H] since we used log gf = —3.93 from iKuruc d {l993). 
However, for our stellar sample we see that the general ap- 
pearance of the [C/Fe] versus [Fe/H] trend is not affected by 
the blend (apart from a somewhat smaller spread in [C/Fe] 
when neglecting it). In what follows we therefore continue 
to use our carbon abundances that we derived with the Fe I 
blend and these are the values we hereafter refer to unless 
otherwise indicated. 

4.5 Error analysis 

To check for systematic effects we compare our derived 
carbon abundances to other studies that also use the 
forbidden [C i] lin e as indicator of the carbon abu n dance 
iGustafsson et ail Il999l : lAndersson fe Edvardssonl Il994t 
iTakeda fe Hondal I2005D In Fig. 03 we plot [C/H] versus 
[Fe/H] for our data set as well as for data from these 
three studies. For each data set we have done a simple 
linear regression between [C/H] and [Fe/H]. These are 
over- plotted in Fig. 03 The regression li nes are given by 
the following equ ations (where G99 =|Gusta fsson et al l 
119991: AE94=lAndersson fe Edvardssonl ~ Il994t 

TH05= lTakeda fe Hondall2005l) : 

[C/H] = [Fe/H] ■ 0.83 - 0.00 (Thiswork), (1) 

[C/H] = [Fe/H] -0.87 + 0.01 (Thiswork, noblend),(2) 



[C/H] = [Fe/H] ■ 0.82 + 0.07 (G99), (3) 

[C/H] = [Fc/H] ■ 0.72 + 0.04 (AE94), (4) 

[C/H] = [Fc/H] ■ 0.82 - 0.12 (TH05). (5) 

The tight relation between [C/H] and [Fe/H] that we 
see for our stellar samp le (see Fig. |SJ is also present in the 
iGustafsson et all lll999l) data set (see also their Fig. 3). As 
can be seen the slopes (A[C/H]/A[ Fe/H]) are essentially the 
sam e for our stellar sample a nd the IGustafsson et alJ (Il999t) 
and ITakeda fe Hondal (|2005) samples. There are, however, 
offsets present between these two data sets and ours. Rel- 
ative to our relationship, the Gusta fsson et ail jl999l) re- 
lationship is shifted 0.07 dex upwards (to higher [C/H] val- 
ues) while the lTakeda fe Hondal J2005I) relationship is shifted 
0.12 dex down wards (to lower [C/H] va l ues). The relation- 
ship from the lAndersson fe Edvardssonl jl994l) data set 4 is 
shifted upwards by 0.04 dex relative to our relationship and 
diffe rs from the other three in that it s hows a shallower slope. 
The lAndersson fe Edvardssonl Jl994l) and ITakeda fe Hondal 
fcooa)" data sets have a conside rably larger spread. In the 
lAndersson fe Edvardssonl ll 19941) case the large scatter is 
due to that their spectra was severely affected by fringing 
and that many of their spectra had an insufficient S/N to 
enable good abundance determination from the extremely 
weak [C i] line. The r eason for the large scatter in the 
ITakeda fe Hondal ll2005l) data is due to that their spectra in 
many cases were of low quality (Honda, private comm.) and 
in many weak-line cases they decided to adopt the values 
from the permitted C I lines (average of AA505. 2/538.0 nm 
lines) as their final carbon abundance instead (note: Fig. |H| 
and Eq. |S]only include data from the forbidden line). 

It thus appear that it is safe to say that there is 
a systematic shif t betw een our abundances and those by 
IGustafsson et all 1^)99) since both data sets are individ- 
ually well defined. Whether this shift is due to the carbon 
abundances (vertical shift in Fig.|HJ or to the Fe abundances 
(horizontal shift in Fig.|SJ, a combination of both, or due to 
that we included the Fel blend and they did not, will now 
be investigated by a comparison of individual abundances 
for stars that we have in common. 

The left hand panel in Fig. [5] shows a detailed 
comparison for 1 3 sta rs that we have in common with 
IGustafsson et al.l ||1999D . The mean difference in abundances 
and stellar parameters are (our values minus theirs, the 
uncertainties represent the one-sigma standard deviations): 
A [C/H] = -0.02 ± 0.07 dex, A[C/H] N o blend = -0.03 ± 
0.07 dex, A [Fe/H] = +0.03 ± 0.06 dex, AT cft = +12 + 66K, 
Alogg = —0.00 ± 0.10 dex. Overall the mean differences 
between our study and IGustafsson et all (|l999() are small 
for all these parameters. It should be noted that for these 
stars our inclusion of the Fel blend barely has any effect 
on A [C/H]. The differences we see in T cS (+12 K) and 
logg (0.00 dex) are small and would only translate into 
A [C/H] » 0.00 and A [Fe/H] » 0.01 (see Table 0. 

In Fig. |5] (right hand panel) we also plot the dif- 
ferences between our abundances and those from other 
works for a total of 31 stars. Comparisons are made 

4 Man y of the [C/H] values in the lAndersson fc Edvardssonl 
ll994l) data set are only upper limit estimates. These were re- 
jected from the analysis presented here. 
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Figure 9. Left hand panel: A detailed comparison of abundances and stellar parameters for stars that we have in common with 
Gustafsson et al. (1999). For the A[C/H] comparison we also include our abundances derived when neglecting the Fel blend (open 
circles). On the "x-axis" we have marked the Hipparcos numbers of the stars. The stars have been sorted accorded to their metallicities, 
[Fe/H] increasing from the left to the right. Right hand panel: Comparison of carbon abundances for stars that we have in common 
with other studies. The figure shows our [C/H] values minus [C/H] values from the other studies and comparisons are made both with 
studies that use the forbidden [C i] line as well as studies that use the permitted C I lines and the molecular CH bands (as indicated in 
the figure). In the case of Andersson & Edvardsson (1994) the arrows mark stars that they only could give an upper estimate on [C/H]. 
In the case of Takeda & Honda (2005) the different symbols mark [C/H] abundances from the forbidden [Cl] line (solid circles) and from 
the permitted Cl lines (open squares). The stars have been sorted accorded to their metallicities, [Fe/H] increasing to the right. On the 
"x-axis" we have marked the Hipparcos numbers of the stars. Note the different scales on the "y-axes" in the left and right hand plots. 
The abundances from the different studies have not been adjusted for differences in the stellar parameters. 



to t he two other studies that use the forbidden [C i ] 
line jAndersson fc Edvardssonll994ilTakeda fc Hondal2005h 
as well as to studies that use the permitted C I lines 
jTomkin et alJl995ilEcuvillon et al.l20q4HTakeda fc Hondal 
20051) and the molecular CH bands iLairdl Il985l) . With 
Andersson fc Edvardsson! <ll994t) we have 12 stars in com- 
mon, and for 4 of those they only give upper limit on [C/H]. 
For the other 8 stars the difference is +0.12 ± 0.12 dex, 
i.e., quite large and substantially larger t han r elative to 
iGustafsson et al] ll99Sh . iTakeda fc Hondal (120051) use both 
the forbidden and the permitted carbon lines and we have 
11 stars in common. The mean difference between our and 
their forbidden [C/H] abundances is +0.04 ± 0.14 dex and 
to th eir permitted [C/H] a bundances it is —0.03 ± 0.07 dex. 
With lTomkin et all (1995) we have 9 stars in common and 
the d ifference is +0.05 ± 0.18 dex. With lEcuvillon et al] 
(2004) we have 13 star s in common and the difference is 
-0.08 ± 0.08 dex. With iLairdl dl985l) we have 11 stars in 
common and the difference is —0.08 ± 0.19 dex 

We would suggest that if on e were to merge our sample 
with the lGustafsson et al] (Il999l) sample one should increase 
their [Fe/H] by 0.03 dex and decrease their [C/H] by 0.02 dex 
(or instead adjust our values accordingly). Even though the 
reason(s) for these shifts are not fully resolved this would at 
least bring our and their s amples onto a common baseline . 
Merging our data with the I Andersson fc Edvardsson! dl994l) 



or the lTakeda fc Hondal ]2005l) samples would introduce a 
lot of extra scatter and make any interpretation of the ob- 
served carbon trends difficult. 

Random errors are partly represented by uncertainties 
in the stellar atmosphere parameters. In Table we list 
how typical uncertainties in effective temperature, metal- 
licity, and surface gravity (AT cff = +70 K, A [Fe/H] = +0.1, 
and Alogg = +0.1, respectively) effects the derived abun- 
dances for three stars (see also discussions in lBensbv et al] 
l2003tl2004l) . 

Random errors may also be due to uncertainties, or de- 
fects, in the observed spectra. Here we paid special attention 
to features in the spectra that could be left-overs from the 
reduction procedure and a possible insufficient removal of 
the fringing pattern. For seven stars (listed in Sect. 14. 31 we 
did not determine any carbon abundances due to that the 
[C i] line was severely deformed. For the remaining stars we 
could not distinguish any possible fringing residuals from 
the ever-present random scatter (noise). If the initial fring- 
ing pattern effects the final carbon abundances (by an ineffi- 
cient removal) we would probably see a substantially larger 
scatter in the carbon trends than what we see in our oxy- 
gen trends for which we used the same spectrograph and 
the same an alysis method but for the forbidden [O i] line 
at 630.0 nm feensbv et al]l2004l) . Although a larger scatter 
indeed seem to be present for the carbon trends than for 
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Table 3. Estimates of the effects on the derived abundances due 
to errors in the atmosphere parameters. <r C ont is the estimated 
error due to erroneous placement of the continuum (see text) and 
(Ttot is the total random error that were calculated assuming the 
individual errors to be uncorrelated. The values for O and Fe were 
taken from Bensby et al. (2004b). 



HIP 88622, (-0.46, 5720, 4.35) 



AT cff = +70 K 
Alogg = +0.1 
Aft = +0.15 km s" 
A[Fe/H] = +0.1 



+0.06 
-0.01 
-0.02 
+0.01 

0.06 



-0.02 
-0.04 
0.00 
-0.03 

0.05 



-0.04 
^0.05 
^0.02 
h0.02 

0.07 



-0.01 
-0.03 
0.00 
-0.01 

0.03 



-0.07 
^0.04 
^0.02 
0.00 

0.08 



-0.03 
-0.01 
0.00 
-0.02 

0.04 
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fO.Ol 
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0.00 


+0.03 


0.00 
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0.00 


A[Fe/H] 


= +0.1 
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+0.03 
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-0.02 
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0.06 


0.06 


0.08 
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= +0.1 
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0.00 
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0.00 
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+0.01 


+0.01 
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0.06 
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0.07 
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our oxygen trends we find that the [C/H] vs [Fe/H] trend is 
quite well-defined (see Fig. |SJ and that effects due to resid- 
uals from the fringing pattern have been reduced to a level 
where they most likely do not effect the general appearance 
of the abundance trends. 

Random errors may also arise from errors in the actual 
fitting of the synthetic spectra to the observed ones. A ma- 
jor concern here is the location of the continuum. Due to the 
weakness of the [C i] line a small change in the location of 
the continuum will have an effect on the carbon abundance. 
This problem will vary with the quality of the spectra as 
well as the strengths of the [C i] line and surrounding lines 
making a true estimate of the uncertainties in the placement 
of the continuum level difficult. By examining our spectra 
and varying the continuum level we estimate an average un- 
ce rtainty of 0.06 dexin [C /H] . Since the spectra we used 
in iBensbv et al.l (120041) did not suffer from fringing effects 
and generally were of higher quality we estimate the error 
in [O/H] to be half of that in [C/H], i.e. 0.03 dex. 

Including a random error of 0.06 dex in [C/H] and 
0.03 dex in [O/H] due to uncertainties in the continuum level 
we estimate that the total random errors in our final [C/H], 
[C/Fe], [O/H], and [C/O] are 0.07, 0.10, 0.07, and 0.07 dex, 
respectively (see Table [3J. 



Table 4. Kinematic properties of the stellar populations that 
were used when calculating probabilities to select thin and thick 
disk stars. The properties of the Hercules stream have been 
adopted from Famaey et al. (2005) and the properties of the other 
populations are the same as in Bensby et al. (2005) except for the 
thin disk that has a lower number density (X) due to the now 
included Hercules stream. All velocities and velocity dispersions 
are given in kms . 





X 


<TTJ 


(TV 






Mag 


Wlag 


Thin disk 


0.84 


35 


20 


16 





-15 





Thick disk 


0.10 


67 


38 


35 





-46 





Hercules 


0.06 


26 


9 


17 


-40 


-50 


-7 


Halo 


0.015 


160 


90 


90 





-220 






5 RESULTS AND DISCUSSION 

Before we discuss our carbon results and where carbon might 
be produced we first, briefly, recap some of the most impor- 
tant features of our sample selection, i.e., the kinematic def- 
inition of the two samples. We then turn to a description of 
the abundance trends found and their implications for where 
carbon can have formed. 



5.1 The kinematic definition of our samples 

The stars in our sample were selected, based on their kine- 
matics, to either be typical representatives of the thin or the 
thick Gal actic disk. Th e kinematic selection is further de- 
scribed in IBensbv et alJ ll2003l . 120051) including an extended 
discussion of how e.g. the local normalisation of the thick 
disk influences the categorisation of an individual star. 

In Fig. HOfa we show a so called Toomre diagram for our 
stars. The distribution of the stars in our two sa mples in this 
plot is essentially the same as lFuhrmanrJ l|2004h finds for his 
thin and thick disk samples. Figure 110b shows the VFlsr 
velocity (which is proportional to the maximum distance 
below/above the Galactic plane, Z m ax, a star can reach) 
as a function of [Fe/H]. We see that even at [Fe/H] = 
or higher there are stars with high Wlsr- The number of 
such stars appear to be higher than one would expect from 
orbit al heating by e.g. molecular clouds llHanninen fe Flvnnl 
2002). One should though remember that our sample is far 
from complete and we certainly are subject to various biases. 
However, our point is that these stars do exist, are readily 
found in any large catalogue (hence not rare objects) and 
they imply that the thick disk indeed extends to at least 
solar metallicities. 

Fi nally, it is also worthwh ile to note, as pointed out 
by e.g. iNordstrom et all l)2004h . that there is a lot of kine- 
matic structure when studying large samples of nearby stars. 
One such structure, or group of stars, of particular interest 
here is the Hercules stream which has kinem atic properties 
very s imilar to some of the thick disk stars. iFamaev et alJ 
(2005) studied a sample of ~ 6700 nearby K and M giants 
and were able to identify several such structures (includ- 
ing the Hercules stream) for which they determined kine- 
matic properties as well as local number densities. From 
their results we see that it is especially Hercules stream 
stars with high C/lsr velocities in the direction away from 
the Galactic centre that could be erroneously classified as 
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Figure 10. a) Toomrc diagram; dotted lines indicate constant peculiar space velocities, i> pe c = (^lsr ^I?SR ^lsr) 1 ^ 2 ' ' n s * e P s °f 50 
km s . b) Wlsr versus [Fe/H], The plots include our full thin disk (empty circles) and thick disk (filled circles) samples from Bensby 
et al. (2003, 2005). Those stars that are included in this study are marked by crosses. 



thick disk stars. Using a co mpilation of chemical abundances 
ISoubiran fc Girardl ll2005l) find that the abundance trends 
for the Hercules stream mostly follow those in the thin disk. 

We have checked for possible contamination from Her- 
cules stream in our samples. This was done by recalculating 
the probabilities that we use d when selecting our samples 
(see lBensbv et al.ll2003t 120051) . and now including the prop- 
erties for the H ercules stream (see Table ffl. Using the same 
equations as in iBensbv et ail ll2003t 12003) (but now modi- 
fied to also include C/i ag and Wi ag ) we find that none of our 
stars are more likely to belong to the Hercules stream than 
to the thin or thick disks. 

For the discussion of our carbon results from this paper 
we will assume that our stellar samples are representative for 
the thin and thick disks, or the stellar components that can 
be kinematically associated with the two disks disks, and 
that the thick disk does extend up and until solar metallic- 
ities. 



5.2 Observed trends 

5.2.1 [C/Fe] versus [Fe/H] 

In Fig. Illfa we show the trend of [C/Fe] versus [Fe/H]. 
The thin and thick disk [C/Fe] trends are fully merged and 
at sub-solar metallicities the [C/Fe] values all fall within 
< [C/Fe] < 0.2 with no particular slope. At super- 
solar metallicities the [C/Fe] values decrease with increasing 
[Fe/H] (independent of if the Fel blend is included or not, 
compare Figs. [TJl and e). 

The flat [C/Fe] trend that we see at [Fe/H] < is 
what generally is found, within the uncertainties, in th e 
literature (see, e . g., co mpilation by iGavilan et alJ |2005|) . 
iGustafsson et all <ll999l) . however, find an increase in [C/Fe] 
for decreasing [Fe/H]. This is somewhat discomforting since 
both our and their carbon abundances are based on the for- 
bidden line at 872.7 nm. Ho wever, at a closer l ook a t the 
[C/Fe] versus [Fe/H] trend in Gustafs son et al.l |l999) it is 
possible that ours and theirs trends are not that different. 
That they see a slope seem to be mainly based on a few (~ 5) 



stars with [C/Fe] > 0.2 at [Fe/H] < -0.6 (see their Fig. 4). 
So, if those stars are disregarded, our and their [C/Fe] versus 
[Fe/H] trends are, within the uncerta inties, similar. Good 
agreement to IGustafsson et all 1^)99) is further strength- 
ened from the detailed comparison in Sect. 14.51 



5.2.2 [C/O] versus [O/H] 

Fig. Illb shows the [C/O] versus [O/H] trend for our stel- 
lar sample. The thin and thick disks show trends that are 
clearly separated. The thin disk shows a shallow increase in 
[C/O] with [O/H] whilst the thick disk first has a flat [C/O] 
trend that increases sharply at [O/H] = 0. The great resem- 
blance with the observed [Fe/O] versus [O/H] trends shown 
in Fig. I11H could indicate that C and Fe indeed originate 
from objects that evolve on similar time scales. As our stellar 
sample spans a rather limited range in [O/H] it is difficult to 
use this sample only to say how the [C/O] trend for the thick 
disk would have evolved from very low oxygen abundances 
(i.e. [O/H] < -0.4). In Fig.[Hlwe therefore show our [C/O] 
versus [O/H] trend together with the stell ar sample (con- 
sisting of stars with halo kinematics) from lAkerman et alJ 
(2004). They based their carbon and oxygen abundances on 
permitted C I and O I lines. As they assumed that the NLTE 
corrections for the C and O abundances for these lines are 
of similar size they did not apply a ny NLTE corrections . 
However, recent NLTE calculations bv lFabbian et all <l2005l) 
have shown that these "permitted" [C/H] values should be 
decreased by as much as ~ 0.4 dex at [O/H] ~ —2.3 and less 
at higher [O/H]. The NLTE corrections for [O/H] are slightly 
less (—0.1 to —0.3 dex, depending on stellar parameters) re- 
sulting in the [C/O] b eing overestimated by ~ 0.2 dex at 
the lowest [O/ H] in the lAkerman et all d2004l) data (see also 
lAsplundll20 05). As a consequence the uptur n seen in [C/O] 
with d ecreasing [O/H] for [O/H] < — 1 in the lAkerman et ail 
(2004) data is overestimated - their trend should be natter. 
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Figure 11. a) and b) show the carbon and oxygen abundances, respectively, relative to iron, c) and d) show the carbon and iron 
abundances, respectively, relative to oxygen. All carbon abundances are from this work and the oxygen and iron abundances have been 
taken from Bensby et al. (2003, 2004, 2005), see also TablelTl Thin and thick disk stars are marked by open and filled circles, respectively. 
The few thick disk stars from Nissen et al. (2002) included in the oxygen plots are marked by filled triangles. 



5.3 Where is carbon made? 

The carbon yields for various objects differ significantly. The 
yields are sensitive to a number of factors such as: stel- 
lar winds, treatment of convection, and the 12 C(a,7) 16 
rate. Earlier models of massive sta rs found them to give 
high yields of C (e.g.. iMaederl 119921) whilst these were ad- 
justed downwards once rotation was taken into account 
jMevnet fc Maederl EoO^) . Several studies have used these 
yields to model galactic chemical evolution and compare 
it with observed abundance trends (both in local stars 
as well as in Galactic and extra-galactic Hn regions). 
iGustafsson et all ill 9991) gave an exhaustive list and discus- 
sion of possible sites for carbon. These include: supernovae, 
novae, Wolf-Rayet stars, low and intermediate mass stars in 
the planetary nebula phase or by super-winds at the end or 
the red-giant phase. Through an analysis of Galactic dwarf 
stars (a mixture of thin disk stars, metal-poor thick disk 
stars, as well as stars with intermediate kinematics from 
Edvardsson et al. 1993) and data from dwarf galaxies they 
arrived at the conclusion that the main contribution of car- 
bon to the galactic chemical evolution comes from super- 
winds of metal-rich massive stars, that low mass stars were 
ruled out, and th at the role of the in termediate mass stars 
remained unclear. iHenrv et alJ (120001) modelled a large set of 
Galactic and extra-galactic H II regions and reached a similar 
conclusion. Other studies find that the low and intermedi- 



ate mass stars are the essential contributors to the chemical 
enrichment of the Galaxy. By adopting the new yields by 
iMevnet fc Maederl ]2002MChiappini et al] {2003b) conclude 
that massive stars do contribute to the carbon enrichment 
but that the main contribution at higher metallicities must 
be due to l ow and intermediate ma ss stars. Using the same 
set of vields lchiappini et al] (l2003aTl are also able to explain 
the abundance ratio s in other galaxies , such as irregulars. 
The recent study bv lCarigi et al] (I2005T) arrives at a similar 
complex explanation for the observed abundance patterns 
where metallicity dependent yields are necessary both for 
the massive stars (increasing yields with metallicity) as well 
as for low and intermediate mass stars (decreasing yields 
with increasing metallicity). In their best-fitting model the 
massive stars dominate the C enrichment at early times and 
the low mass and intermediate mass starting to contribute 
later. The contributi ons at later t i mes a re roughly equal for 
the two s ources. A sJCar igi et alJ (2005) had to invoke the 
yields by IMaederl dl992D ( that" no longer should b e used 
as they are superseded bv lMevnet fc Maederl 120021) in or- 
der to reproduce the solar [C/O] ratio the importance they 
attribute to massive stars at higher metallicities could be 
overestimated. 

Our [C/Fe] vs [Fe/H] trend in the thick disk is flat 
whilst our [O/Fe] vs [Fe/H] trend (see Figs. II lb and b, re- 
spectively) shows a clear break and subsequent downturn at 
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Figure 12. [C/O] versus [O/H], Our thin and thick disk stars are marked by open and filled circles, respectively. Halo stars from 
Akcrman et al. (2004) are marked by crosses. Note that the Akerman et al. (2004) [C/O] values are overestimated by ~ 0.2dex at the 
lowest [O/H] due to neglect of NLTE effects. The effect decreases as [O/H] increases. 



[Fe/H] = -0.5 (this corresponds to [O/H] ^ 0). This type of 
dow nturn is normally interpreted a s the onset of SN la (see, 
e.g.. lTinslevlll97fll : lMatteuccill200ll) and since oxygen is only 
produced in SN II the increased Fe production by SN la re- 
sults in a downward trend as [Fe/H] increases. For [C/Fe] vs 
[Fe/H] we see no such trend. Hence, we may infer from the C 
abundances in our thick disk sample that enrichment from 
carbon happens on the same time scale as the enrichment 
from SN la. Hence the increase in Fe production is matched 
by enrichment of C. 

The fact that our [C/O] -[O/H] and [Fe/O] - [O/H] 
trends are indeed very similar further strengthens the im- 
portance of low and intermediate mass stars as contributors 
to the carbon enrichment at higher met allicities. This is also 
supported bv IChiappini et al.H200 3b'l who predict, if low 
and intermediate mass stars are important, that there first 
would be a flat [C/O] plateau for the thick disk that then (at 
roughly solar [O/H]) should sharply increase, to be followed 
by a shallow thin disk [C/O] trend at higher [C/O] ratios. 
And this is j ust what our data show . 

Further, lAkerman et al.l i2004t) note that the ampli- 
tude of the rise that they see in the [C/O] ratio around 
[O/H] = -1 in their Milky Way halo stars (see Fig. is 
well matched by similar trends for H II regions in nearby spi- 
ral and irregular galaxies. This they interpret as favourable 
for their findings that massive stars are the main contribu- 
tor to the chemical enrichment of carbon as it would be very 
unlikely that the abundance trends in galaxies that have ex- 
perienced different types of chemical histories should be so 
similar if carbon was mainly produced in sites where a sub- 
stantial time-delay had to be invoked, e.g. low-mass stars. 
Given the uncertainti es relating to NLTE corrections of the 
AJ«rman^t^Ll (|2004jWlata (see discussion above) and that 
Chiappini et al (2003a) are able to explain the abundance 
ratios in other galaxies as well (including irregulars), using 
the same set of yields that were used for the Milky Way, 
this indicate s that the role o f mass ive stars might be over- 
estimated bv lAkerman et ail i2004l) . 

If we compare our C abundances with Y the trend 
of [C/Y] versus [Y/H] shows a flat trend for [Y/H] < 
which turns to a steadily declining trend for [Y/H] > (see 
Fig. 113b ) which is very similar to what we see for [C/Fe] 
versus [Fe/H] (see Fig. Illfa). [C/Y] vs [Fe/H] is somewhat 



different, being, essentially, flat for both the thin and the 
thick disk samples at all metallicities (see Fig. 113b ). Approx- 
imately 70% of Y is produced in the s-process in AGB stars 
and ~ 30% comes from massive stars (~ 10% from a weak 
"secondary" s-pr ocess component and ~ 20% from a pri- 
mary component iTravaglio et~ai]|2004l) . The most straight- 
forward interpretation of the lack of trends would be that 
the two elements, C and Y, are made in objects that en- 
rich the interstellar medium on the same time scale. Most 
naively we would assume that their major components are 
indeed made in the same objects, namely low and interme- 
diate mass stars in the AGB phase. However, if we inspect 
the trend of [C/Y] versus [Y/H] in addition to the flat trend 
below sub-solar [Y/H] we see a declining trend of [C/Y] for 
[Y/H]> indicating that as the Y increases there is a de- 
crease in C production. This could be due to a metallicity 
effect in the production of s-process elements that favours 
light s-process elements (such as Y) over heavy s-process el- 
ements (such as Ba ) at high metallicities ( Buss o et al 1l200ll : 
iBensbv et aT]l2005l) . 

So, we have shown that the carbon production must bal- 
ance the the Fe production by SNIa and the Y production 
by AGB stars. With this observable in mind it appears that 
either the massive stars must have (very) metallicity depen- 
dent yields that are finely tuned to the [Fe/H] timescale of 
SNIa as observed in the Milky Way thick disk or that less 
massive stars are significant contributors too. 

In summary we thus find that the carbon enrichment 
at low metallicities (i.e. the Galactic halo and metal-poor 
thick disk) is due to massive stars but as more and more 
low and intermediate mass stars start to contribute to the 
general enrichment they also more and more dominate the 
carbon enrichment of the interstellar medium as it is seen in 
the thin disk and the metal-rich thick disk. 



6 SUMMARY 

We present carbon abundances in 51 F and G dwarf stars 
(16 thick disk stars and 35 thin disk stars)in the solar neigh- 
bourhood. The analysis is based on the forbidden [C 1] line 
at 872.7 nm which is an abundance indicator that is insensi- 
tive to errors in the stellar atmosphere parameters. Combin- 
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From our analysis of the forbidden [C i] line at 872.7 nm 
to determine carbon abundances in solar-type stars we also 
found a few details that we think are noteworthy: 1) By 
examining how the neglect of the blending Fel line effects 
the derived carbon abundances from the forbidden [C i] line 
we find that it is mainly for stars with effective temper- 
atures less than ~ 5700 K that differences are large (see 
Fig- Ek); 2) Currently, the only available source for the os- 
cillator streng th of the blending Fe I line is a theoretical one 
iKurucz)ll993h . Therefore, it is possible that the strength of 
the Fel could under- or overestimated. By working strictly 
relative to the Sun the effects of an erroneous treatment 
of the blend will be reduced. It would, however, be very 
valuable if a log gf -value for the Fel line could be measured 
in the laboratory; 3) For the Sun we find a carbon abun- 
dance of £q(C) = 8.41 when including the Fel blend and 
60(C) = 8.44 when neglecting it; 4) The general appearance 
of our carbon trends in the two disks is similar whether we 
include the Fe I blend in the abundance analysis or not (see 
Figs. |3l and e). 



Figure 13. a) [C/Y] versus [Y/H] and b) [C/Y] versus [Fe/H], 
Yttrium abundances are taken from Bensby et al. (2005). Thin 
and thick disk stars are marked by open and filled circles, respec- 
tively. 



ing these data with our previously published oxygen abun- 
dances, based on t he forbidden [0 1] line at 630.0 nm (see 
iBensbv et al"ll2004h . we can form very robust [C/O] ratios 
that we then used to investigate the origin of carbon and the 
chemical evolution of the Galactic thin and thick disks. Our 
findings concerning carbon abundance trends in the Galac- 
tic thin and thick disks include: 

* At sub-solar [Fe/H] we find that the [C/Fe] versus [Fe/H] 
trends for the thin and thick disks are fully merged and 
essentially flat. For the thin disk, that extends to higher 
metallicities, [C/Fe] is flat until [Fe/H] = when a slow de- 
cline starts (see Fig. Illb l. 

~k Our abundance trends indicate that the sources that con- 
tribute to the carbon enrichment of the interstellar medium 
do so on the same time-scale as those that produce most 
of the iron, i.e. SNIa. The production of C and Y (com- 
ing mainly from AGB stars) also seem to work on the same 
time-scale. 

it In light of our own as wel l as other studies in the l iteratu re 
(notably iGustafsson et all Il999l: IChiappini et al.l l2003aHbt 
lAkerman et all 120041: ICarigi et alJ 12003) we feel that the 
source(s) of carbon is not yet settled but that there is grow- 
ing evidence that a complicated, and finely tuned, set of 
objects contribute to the enrichment of carbon in galaxies. 
As discussed, based on our own results only we would con- 
clude that the main source for carbon in the Galaxy is low 
and intermediate mass stars. However, for the Milky Way 
galaxy it appears that massive stars played a significant role 
for the carbon enrichment at low metallicities (i.e. halo and 
metal-poor thick disk) whereas low and intermediate mass 
stars dominate more and more at higher metallicities, i.e. 
that they have been the major contributors to the carbon 
enrichment in the thin disk and the metal-rich thick disk. 
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